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Abstract. Fleming's inequality is generalized to the decay function of mixed 
states. We show that for any symmetric hamiltonian h and for any density 
operator p on a finite dimensional Hilbert space with the orthogonal projection 
n onto the range of p there holds the estimate Tr(ne-"'*pe"") > cos'^ {{Ah) pt) 
for all real t with {Ah)p \t\ < ■7t/2. We show that equality either holds for all 
t € M or it does not hold for a single t with < {Ah)p \t\ < -k 12. All the density 
operators saturating the bound for all t G R, i.e. the mixed intelligent states, are 
determined. 



PACS numbers: 03.65.-w 

1. Introduction 

Two states p\ and pi of a quantum system can be discriminated on the basis of a 
single measurement outcome if there exists an observable A such that the probability 
measures which are generated by p\ and pi on the spectrum of A have disjoint sup- 
ports. In particular if a state p evolves under a Hamiltonian H into the state pt it 
may be desirable to determine and perhaps to minimize a time t > when the evolved 
state Pt can be discriminated from the initial state p by a single measurement. A more 
realistic goal is to distinguish pt from p by performing single measurements on "few" 
ensemble members only. If one chooses as observable A an orthogonal projection 11 
with Tr (Hp) = 1, then this can be done if Tr (Xlpt) is close to since this means that 
it is very unlikely to find the property 11 in the state pt, while it is certain in the state p. 

Under a somewhat broader perspective the quantity Tr(n/9t) is commonly used 
in order to formalize the intuitive picture of the decay of a property. p[ Since in many 
cases the survival probability Tr (Ilpf) of the property 11 cannot be computed explic- 
itly, there arises the quest for estimates of the decay- function Pp : 1 1— > Tr (Ilpt) . 

An important such estimate for Pp in the case of a pure state p and in case of 
the property 11 = p is due to Mandelstam and Tamm This estimate has been 
rediscovered by a different reasoning almost 30 years later by Fleming. 3 Since then 
it is called Fleming's bound. It says that for any pure state p with a finite energy 
uncertainty (AiJ)^ there holds 

Pp (t) > cos^ ^ — --^ for aU t with < 7r/2. (1) 
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From the estimate ([T|) a lower bound to any positive t such that Tr (Ilpt) = e is 



The special case e = leads to the inequality 



for the smallest time t > with Tr (Upt) — 0. This time, if existent, is called orthog- 
onalization j4] or passage time [5]. Clearly it would also be useful to have an upper 
bound for Pp, from which the existence of an orthogonalization time could be inferred. 
Polynomial upper bounds have been given by Andrews jOJ, which, however, are strictly 
positive. Therefore they do not yield an upper bound to an orthogonalization time. 

A simple geometric meaning of Fleming's bound became clear through the time- 
energy uncertainty relation of Aharonov and Anandan [7j- First, 2t{AH)p/h equals 
the arc length of the curve X i-^ p\ with < A < t in the projective space V (H) of 
one dimensional subspaces of H. Second, 2 arccos ^/Pp(J) equals the geodesic distance 
between p and pt in V (7i) . Here the Riemannian geometry is defined by the Fubini- 
Study metric of V {Ti) . Thus, as has been pointed out by Brody [5], Fleming's bound 
([T]) is equivalent to the fact that the length of a curve in V {H) is not less than the 
geodesic distance between its initial and end point. 

In Ref. [S] for given Hamiltonian H all pure states p with an orthogonalization 
time equal to the lower bound Trh/2 {AH)^ of equation ^ have been identified. I.e., 
for such states there holds {AH)pt — hn/2 for the smallest t > Q with Pp [t) ~ 0. 
These states are called "intelligent states" as they saturate the Aharonov-Anandan 
uncertainty relation. A pure state p is found to be intelligent if and only if there 
exist two eigenvectors 0i ,4)2 oi H corresponding to different eigenvalues and with 
W'PiW = \\'t'2\\ such that p equals the orthogonal projection onto the one dimensional 
subspace <C ■ {^i + ■ [8] 

In [5] a lower bound for the smallest t > Q with Pp (t) = is given for a special 
type of a mixed states: The density operator p is assumed to be a mixture of mutually 
orthogonal intelligent pure states. In Ref. [1] a generalization of the orthogonaliza- 
tion time to mixed states has been addressed too. In this work, however, the fidelity 
Tr \/ \fppt\fp is used as a measure of the degree of decay. Clearly, for mixed states 
the fidelity does not coincide with Tr (Ilpt) . So neither of the two works [5], and [4] 
presents a generalization of Fleming's bound ([T]) to the case of an arbitrary mixed state. 

In this paper we generalize Fleming's bound to an arbitrary mixed state p. We 
consider the decay function Pp (t) — Tr (Ilpt) , where 11 is chosen to be the orthogonal 
projection onto the range of p and we confine ourselves to finite dimensional Hilbert 
spaces. We first extend Fleming's bound to Pp. We then sharpen the bound by proving 
that only one of the two cases 




arccos ^fe < t. 



m)p 




I Clearly this does not imply that there exists any t at all such that Pp {t) = e holds. 
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(ii) Pp (t) = cos^ t/hj for alH e M 

is realized. Then we identify the set of ah density operators which saturate 
Fleming's bound. In order to have the paper reasonably selfcontained we have included 
a treatment of some closely related well known results on pure state decay. In this 
way it also becomes more visible which structures remain unchanged when going from 
pure states to mixed ones. The pure state decay function Pp is denoted as P^ when 

2. Pure state decay 

Let H be a finite dimensional Hilbert space. The scalar product of two vectors (j),tp gH 
is denoted by (</>, V') • Let the dynamics of H be given in terms of a symmetric linear 
operator h : TC ^ H hy (pt = cxp (—iht) (/) for f e R and (f) E H. The survival amplitude 
: R ^ C is defined for e H with = 1 through (t) = {(f), cj)t) and accordingly 
the survival probability of as a function of t is given hy P^ — \A^\ : R ^ R>o- 
From the Cauchy-Schwarz inequality we have < 1. The nonnegative number [t) 
is the probability that the pure state (j)t (</>t, •) passes a preparatory filter for the state 
(j) {4>, •) . Due to 



A^ {-t) = A^ {t) 

P^ is an even function. Since (po = (p we have A^ (0) = 1 = P^ (0) . 

The expectation value of h in the state </> (</>, •) is denoted by {h)^ = (0, h(p) and 
its variance reads 

{Ah)l = (h')^-{h)l. 

4> is an eigenvector of h if and only if (A/i)^ = 0. Thus for (A/i)^ = the function P^ 
is constant, i.e. {t) = 1 holds for all t. For (A/i)^ > 0, however, P^ is not constant 
since for f ^ 

l-it{h),-lt^{h^)^ + il^t^h^)^ + 0{t')' 



P<t>it) 



= l-{/^h)lt^ + 0{i^). 
Thus P^ has a strict local maximum at t = if and only if (A/i)^ > 0. 

Due to the spectral theorem there exist (unique) nonzero pairwise orthogonal 
vectors <pi, . . .(f)„ with h(j)a — i^a4>a and wi < . . . < a;„ such that 

cj)t = e-'"i Vi + • • • + e-'""Vn 

for all t. Then A^ {t) = YZ=i Aae"''^"* with A„ = ||0af > follows. For P^ {t) one 
obtains 

n n 
a, ,8=1 Q:,/3=1 
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Thus both and are the restriction of an entire function to the real Une. In 
particular and are C°° functions. 

It has been shown by Mandelstam and Tamm [5], and along a different strategy 
by Fleming in [3J that for all t with (Ah). \t\ < 7r/2 there holds 



>cos2((A/i)^t) 



The original proof of Mandelstam and Tamm [2] has been elaborated by Schulmann 
in [9]. A new proof has been given recently by Kosihski and Zych. [TO] 

We shall now prove the following somewhat stronger result implicitly contained 
in [4] and [5]. 

Proposition 1 Let (f> E H with \\<j)\\ = 1 and (Aft,)^ > 0. Then exactly one of the 
alternatives (i) and (ii) holds. 

(i) P0 {t) > cos^ (^(Ah)^ t'j for allt€R with < (Ah)^ \t\ < 7r/2 

(u) {t) = cos^ {{Ah)^ t^ for allteR 

Alternative (ii) holds if and only if there exist two vectors (/>i,^2 G with 
h(j)i = uJi(j)i, LUi < LU2, ||0i|| =1/2 such that (j) = + 02- 

Proof. LetU = (j) {(f), ■) . Then holds P^ (t) = (0, e^^'He-''^'^) = (H)^^ . From this it 
follows that 

(t) = i (0, e"'* [h, n] e-"'*^) =i{[h,U])^^. 

Using the uncertainty relation for the pair {h, II) we thus obtain for P^ (t) = -^P^, (t) 
the estimate 



\p^{t)\ = {[h,n])^^ < 2 (A/i)^ (An), 



From (An)^^ = (H^) - {n)l = (H)^^ - {n)l = (n). 1 - (n>^J it follows that 



for alH e M 

|p;(t)| <2(A/i)^y^P^(t)(l-P^(t)). (4) 

We first simplify this inequality by introducing the dimensionless time variable 
X — t {Ah}^ and the function w : M [0, 1] with v (x) — P^ (i) . Inequality ([4]) then 
becomes equivalent to 

-2y^v {x){l-v{x)) < v' {x) < 2^Jv (x)(\-v(x)) for aU x G M. 

In order to make use of the differential inequality 



v' (x) (5) 
we first discuss the (autonomous) differential equation 



y' = f (x, y) with / : M X (0, 1) ^ R, / [x, y) = -2^y[l-y). (6) 

The function yo '■ (0:i'/2) (0,1) with yo{x) ~ cos^ x is a solution of this 
differential equation since for all x G (0,7r/2) 



y'o i^) = -2 cos (x) sin (x) = -2y/yo (x)a/1 - yo (x) = / (x, yo (x)) . 



Fleming 's bound for the decay of mixed states 



5 



This solution of ^ is of maximal domain since the limits 
lim yo (x) = 1 and lim yo (x) = 




do not belong the admitted range < j/ < 1 of solutions. Other solutions of maximal 
domain are obtained from yg by translation: yc (x) — yo {x — c) for c < x < c + 7r/2. 
By a suitable choice of c the initial value problem j/c {£,) = rj for any 77) g M x (0, 1) 
is solved. Since / obeys the local Lipschitz condition of the uniqueness theorem for 
the solutions of first order differential equations, the set of all solutions to y' — f (x, y) 
with maximal domain is given by {ydc M} . 

The continuous extension g of / to the domain R x [0, 1] leads to the differential 
equation z' = g {x, z) = —l^z (1 — z) which violates the local Lipschitz condition on 
the boundary points (a;, z) with either z = or z = 1. The set of solutions of the 
extended equation with maximal domain is given by {zc |c G M} with 



Thus any function z^ with c > is a solution of the initial value problem z (0) = 1 
with maximal domain. For any such solution z^ with c > holds 

Zo {x) < Zc (x) < 1 

for all x > 0. 

According to a theorem of differential inequalities, quoted in the appendix A, we 
then conclude from ^ and from i> (0) = 1 that for all a; > 

V (x) > Zq (x) . 

Thus V (x) > cos^ X for all x E [0, 7r/2] . This is Fleming's inequality. 

Suppose now that 77 :— v{^) > cos^ ^ for some ^ G (0, tt/2) . With rj = cos^ — c) 
for some c G (0,7r/2) it then follows again from the quoted theorem on differential 
inequalities that v (x) > cos^ {x — c) > cos^ (x) for all x G [£,, 7r/2] . From Fleming's 
inequality we now have the two cases only: 

(i) For any e > there exists a ^ G (0, e) with v{^) > cos^ ^. 

(ii) There exists an e > with v (x) ~ cos^ x for all x G (0, e) . 

In the case (i) we have v (x) > cos^ (x — c) > cos^ (x) for all x G [^,7r/2]. 
Since there exist such ^ arbitrarily close to it follows that v (x) > cos^ (x) for 
all X G (0, 7r/2]. Since v is an even function the inequality extends to all x with 



In case of (ii) the identity theorem of holomorphic functions implies v (x) = 
cos^ (x) for all X G M since v is the restriction of an entire function to the real line. 
Thus we have derived the alternatives (i) and (ii) as being exhaustive. 




1 for X < c 

cos^ (x — c) for c < X < c + 7r/2 
for X > c + tt/2 



\x\ G {0,tt/2 
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Suppose now that alternative (ii) holds. From the spectral decomposition ([3]) of 
Pa we extract the constant term and the one with the highest frequency according to 



a,/3=l 
Q>/3 



^ + 2A„Ai cos [{ujn - wi) t] + 2 ^ AqA/j cos [(wq - ojp) t] . 



a, 13=1 
Q>/3,(a,/3)^(n,l) 



The assumption {t) — cos^ (^(Ah)^ t^ = ^ ^1 + cos ^2 (Ah)^ t^^ now implies, due 
to AqA^ > for all a,f3, that the index set of the last sum is empty. Thus we have 



n = 2 and 



A? + A^ = i, 2A1A2 = ^, - wi = 2 {Ah)^ 



The first two equations imply Ai = A2 = 1/2. From this it follows that the third 
condition a;2 — = 2 (Ah)^ holds, since 

(A/l)^ = XiUjf + \2'jj\ - (AlCJi + A2W2)^ 

^ \ (^1+^2) - \ {^l + ojlf 

= ^(wi -W2)^ 

Thus we have derived from alternative (ii) that is a linear combination of 
just two eigenvectors of h with spectral components of equal norm. The inverse 
conclusion that alternative (ii) follows from </> = </>! +02 with h(j)i = LiJi(j)i,<jj2 > wi and 
||0i||^ = Ai = 1/2 is obvious from 

P<P (t) = A? + A2 + 2A1A2 cos [(cj2 - wi) i] 

= i (1 + cos [(W2 - OJl)t]) = cos^ ((A/i)^ . 



3. Mixed state decay 

Let p : TC ^ TC he a density operator on the finite dimensional Hilbert space H, 
i.e. p is linear with p > and Tr (p) = 1. Due to the spectral theorem there exist 
mutually orthogonal vectors ipi, . . .ipn with ||V'fe|| — 1 for all k and there exist numbers 
Ai, . . . A„ S M>o with X]fc=i Afe = 1 such that 

n 

P = ^>^kil'k{ipk,-) ■ (7) 

k=l 

The orthogonal projection Ii : H ^ H onto the range of p is given by 

n 

n = ^Vfc(V'fc,-). 

k=l 

The projection II is the smallest orthogonal projection with Tr (pll) = 1. 
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For an arbitrary orthogonal projection E : H ^ H the nonnegative number 

Tr (pE) is the probabihty that the state p passes a fiher for the property associ- 
ated with E. More generally, the expectation value of a linear symmetric operator 
A -.n^n is given by (A) ^ = Tr {Ap) and its variance reads (AA)^ = (A^)^ - (A)^ . 

The dynamics = exp {—iht) (j) is extended from vectors to density 

operators through p\-^ pt = e~''**/oe'''*. As a generalization of the survival probability 
to mixed states one may consider the function Pp -.M. ^ [0, 1] with 

Pp (t) = Tr (ne-"'Ve"'*) = (e'''*ne-'''*)^ = (H)^^ . 

The number Pp (t) thus gives the probability that the evolved state pt has the property 
n associated with the initial state p. Again t = is an absolute maximum of Pp since 
Pp (0) = 1. Prom this it follows that Pp (0) = since Pp is differentiable. 

Let $1, . . . $g with g > n be an orthonormal basis of H such that h^r = w^^r 
for r = 1, . . . g. Then holds 

9 

Pp {t) = Tr (e'^^'ne-^'V) = e"'*ne-"**p$^> 

r=l 

q 

= ^ e'('^'-'^=)*($„n$,)($„p$^). 

r,5=l 

Thus Pp is a finite linear combination of exponentials and thus of C°° type. 

As in the case of pure states the condition (A/i)^ — implies Pp (t) = 1 for all t. 
This can be seen as follows 

o={Ah)i={h')^-{h)i=(^[h-{h)^y^^ 

Thus we have {h — {h) p)^k = for all k. Therefore all the vectors ?/^fc contributing 
to the spectral decomposition of p are eigenvectors of h (with the same eigenvalue). 
From this then follows the stationarity of p, i.e. Pt = P for all t. While in the case 
of pure states the condition (Ah)^ > implies that P^ is not constant, this is not 
so with mixed states. A counterexample is provided by any p such that 11 commutes 
with h as it is, e.g., the case for p {H) = H, since then 11 = id-n- 

In order to better understand Pp near we first observe 

n 

Pp {t) = TV (ne-^'Ve"^*) = Yl (^fc, e-'^Ve'^'Vfe) 

fe=i 

n n 

k,l=l k,l=l 

n n 

= EAfc|(Vfc,e-^'*Vfe)|'+ J2 ^i\{^k,e-'^%)\\ 

fe=l k,l=l 

kjtl 
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n n 



(8) 



k=l 



The Taylor expansion of Pp at now yields 

n n 

Pp{t) =Y,^k(i-{Ah)ljA+t^ E ^i\{i'k,hi^i)f + o{t^ 



fc=l 



k^l 



fe=i 



k,l=l 
k^l 



From this we infer 



P'JiO) 



(9) 



fc=i 



k^l 



We shall now prove a generalization of Fleming's bound to the survival probability 
of mixed states. 

Proposition 2 Let p : Ti. ^ Ti be a density operator such that {Ah)^ > 0. Then 
exactly one of the alternatives (i) and (ii) holds. 

(i) Pp (t) > cos2 (^{Ah)pt^ for allteR with < (Ah)^ \t\ < tt/2 

(ii) Pp (t) = cos2 (^{Ah)p t^ for allteR 

Alternative (ii) holds if and only if there exist two (different) eigenvalues uji,u!2 
of h such that every vector ipk which appears in the spectral decomposition ^ of p 
has a decomposition ipk — (j)k,i + 4>k.2 with 

h(j)k,l =UJi(f>k,l, h(j)k,2 = UJ2(f>k,2 and {(f>k,e,4'l,ri} = -^h,lSe,r, 

for all k, I G {1, ... 71} and for all £, ?7 G {1, 2} . 

Proof. As in the case of pure states we start from 

it) = ^ (e"'*ne-'''*)^ = i (e"'* [h, U] e""'*)^ ^i{[h, U])^^ . 

The generalized uncertainty relation for the mixed state pt applied to the pair of 
observables (/i,n) reads 

2iAh)^^ (An)^^ > |([/i,n])^^ 

From n-^ = n we obtain (AII)^^ = Pp{t) (1 — Pp (t)) and therefrom the estimate 



Ppit) 



<2{Ah)^ ^Pp{t){\-Pp{t)) 



for all t G 
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The alternatives (i) and (ii) follow from this for i > in exactly the same way as 
in the case of the pure state survival probability P^. Since, however, the mixed state 
survival probability Pp need not be an even function, the case t < needs a separate 
consideration: The case i < is transformed into the case i > by replacing the 
Hamiltonian h through —h. Since the variance of — /i in the state p is the same as that 
of h, the alternatives (i) and (ii) hold for for t < unchanged. 



Suppose now that alternative (ii) holds. Then Pp(t) — cos {{Ah) pt) = 1 — 
{Ah)l + O (<■*) for t 0. Thus -P^' (0) /2 = (A/i)^ holds. From equation ^ 
we then obtain 



{Ah)l = Y,X,{Ah)l^- Xi\{ij,,hi;i)\\ (10) 



k=l k,l=l 
k^l 



Now a general result of probability theory says that the variance of a stochastic 
variable under a mixture of probability measures is greater or equal to the mixture of 
individual variances, or more specifically applied to the present context it says that 



n n 



iAh)l-J2x,{Ah)l^=J2Y.^'^^^{('')^.~(^h.) ^0- (11) 

fe=l fc=l /=fc+l 

The proof of equation pip is given in the appendix. From the equations (jlip . and 
(fTO]) it thus follows that 

71 71 71 2 

o>- 5] \i\{i^kMi)? ^Y. E -W^O ^0- 

kj=l k=ll=k+l 
k^l 

Thus both sides of this equation must vanish and {ipk, hipi) = and {h),^p^ = (^)^, 
follows for all (fc, I) with k I. Furthermore we have 

n 

{^h)l^Y.>^,{Ah)l^. 

k=l 

From jl]) it follows for Pp (t) = cos^ (^{Ah)^ ^) = 5 (l + cos (2 {Ah)^ t^j that 

n 71 

-(l + cos(2(AM^t)) =^AfeP^,(t)+^ A,|(^,,e-"'V/)r.(12) 

fc=l k,l=l 
k^l 

This implies that each of the even functions P^^. is a real linear combination of the 
constant function 1 and cos(2(A/i)pi). Thus we have for alH G M 



P^,, [t) =Ak + Bk cos [2 {Ah)p tj =Ak+Bk- 2Bk sin^ [{Ah)^ t 

= 1 -2PfeSin2 ((A/i)pt) . 

with constants Ak,Bk G M such that P^^ {0) = Ak + Bk = I. From < P^^ (t) < 1 it 
follows that < 2Pfe < 1. 

Thus P^^ obeys for i 

P^, (i)-l-2P, {Ah)y + 0{t*). 
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Taking into account that (V'fe, h4>i) = foi k ^ I the right hand side of equation (fT^ 
obeys 

n n n 

Y,\uP^,{t)+ Ai|(^fc,e-"'*Vi)r = E^fc(l-2^fc(^'^)^') +0 

fc=l k,l=l k=l 

k^l 

Thus we conchide from equation that 

n 

l-(A/.)^t2=^A,(l-2Bfe(A/i)^i2). 

fc=i 

From this it follows that X)fc=i ^k'2,Bk = 1, which in turn implies by means of 
< 2Bk < 1 that 2Bk = 1 for aU k. Thus we have (A/i)^^ = (A/i)^ and 

P^, (0=cos2((A/i)^t) 
for each k. From (fT^ it now follows that 

n 

Ai|(^fc,e-'''Vi)r = 

k,l = l 
k^l 

for all t. For each of the vectors V'fc alternative (ii) of proposition [T] is thus realized. 
From {h)^^ = ^^"^ from (A/i)^^ = (A/i)^ it finally follows that the eigenvalues 

Wfc^e in ^'/'fc.e = ^k,e4'k.e do not depend on fc. 

The inverse statement is obvious by direct computation. ■ 



4. Appendix: Differential inequalities 

Let /, J be two closed real intervals with ry) G / x J and let / : / x J ^ R be 
continuous. Then the following results can be found in either Chapt. I, §9, sects. VI 
and VIII (pp. 73 - 75) of Ref. [11] or in Chapt. II, §8, sects. IX and X (pp. 67 - 69) 
ofRef. [12]. 

Proposition 3 The initial value problem y {S.) = rj of the differential equation 
y' = f (a;, y) has two solutions y<, and y* which both extend to the boundary of I x J 
such that any other solution y of this initial valueproblem obeys {x) < y (x) < y* (x) 
wherever both sides of an inequality are defined]^ 

Proposition 4 Let v : I ^ J and w : I ^ J be functions with 

V {£,) < rj and v (x) < f {x^v (x)) for all x > £^ 

w ^ V '^^'^ ™' (^) ^ / (a;, w (x)) for all x > £, 

then holds v (x) < y* (x) and w (x) > y^ {x) for all x > wherever both sides of an 
inequality are defined. 



§ The solution j/* is called minimal and y* is called maximal. Yet it is also common to call any 
solution of maximal domain a maximal solution. These two notions of maximal solutions thus should 
not be confused. 



Fleming's bound for the decay of mixed states 11 
5. Appendix: Variance and mixing 

Lemma 5 Let p : Ti TL be a density operator on the finite dimensional Hilhert 
space Ti. with its spectral decomposition as given by equation Let h : Ti. s- Ti. be 
linear and symmetric. We abbreviate {h)^,^ by (/i)^, . Then holds 

n n 

(AMp - E ^^^^l + 2 E ^^^^ (^^h - {h)if • 

k=l k,l=l 

Proof. First we observe that 

n n 

{Ah)l = {h')^ {h)l = E A, {h')^ - A.A, {h), {h), 

n n n 

= E A, {Ah)l + E Afe {h)l - J2 >^kXi {h), • 

k=l k=l kj=l 

From the last term we extract the contribution with k — I to obtain for M := 

iAh)i-j:LiXk{Ah)i 

n Ti Ti Tl 

M^Y.^u{h)i-Y,>l{h)i-Y, J2 >^kXi{h),{h),. 

k=l k=l k=ll=l,l^k 

In the second sum of this we replace A^ = A^ ^1 — X^j^^fe ^ij which yields 

n n n n 

n Ti 

E E 

k^i i^i,i^k 

TL TL 

2 

fc=i i^i^i^k 

n n 



E E x,\i({h)i-{h),{h), 

k^i i^i,i^k 

^ n TL 

2E E X,Xi({h)l + {hfi-2{h),{h)^ 

k^l l^l,l^k 

^ n TL 

oE E x,xu{h),-{h)f. 



2 

k=l l=lj=£k 
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